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Abstract. We compute the three loop MS anomalous dimension for the 3-quark operator cor- 
responding to the proton. This requires the treatment of 7 s within dimensional regularization 
as well as evanescent operators generated through the renormalization. We extend the Larin 
scheme for 7 s to a mixing matrix of finite renormalization constants chosen so that chiral sym- 
metry is manifest in four dimensions. We also provide the finite part of the Green's function 
at two loops where the operator is inserted at zero momentum in a quark 3-point function in 
an arbitrary linear covariant gauge in order to assist with the lattice measurement of the same 
quantity. The renormalization of the generalized 3-quark operators in the scheme devised by 
Krankl and Manashov is extended to three loops and the anomalous dimensions for the (^,0), 
(§,0) and (1, 5) spin operators with various chiralities are also given. 
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1 Introduction. 



In quantum field theory baryons are represented by 3-quark operators where the operators are 
chosen so as to have the same discrete and continuous symmetries as the observed states in 
nature guided by the quark model. For instance, protons are built from two up and one down 
quark fields with an overall spin of spin |. With suitable choices of handedness for these three 
fields one can construct a 3-quark operator with the correct J PC values for the proton. While 
in nature these quarks are in general massive and confined, for some theoretical studies of such 
hadronic states one can consider them to be built from chiral or massless fields. This is an 
appropriate approximation in the high energy limit where the masses are small when compared 
to the momentum scale. Thus the leading twist operators will dominate any high energy analysis, 
[II [21 E] . At low energies the internal structure of the baryons and in particular protons, can be 
probed by measuring the structure functions and their moments. For the latter the field theoretic 
quantities of relevance are 3-quark operators involving covariant derivatives, [H El [3j HI El E] . 
Knowledge of the proton structure can assist with understanding the strong interaction in the 
infrared regime and hence how quarks condense or hadronize to form the nucleon states we see in 
nature. In quantum field theory the use of the Quantum Chromodynamics (QCD) Lagrangian 
provides us with a forum to study these different aspects of hadrons. At high energy one can 
apply perturbation theory and calculate order by order in the strong coupling constant which is 
assumed to be small in this region. However, while this provides one with the renormalization 
group evolution of the operators with scale, it is the non-perturbatively measured operator 
matrix element which gives the information relating to the baryon structure. Such matrix 
elements cannot be determined perturbatively for the purpose of extracting physical results. 
Instead they are measured non-perturbatively using lattice gauge theory where the spacetime 
is discretized. Although this involves the use of large computers to handle the huge numerical 
calculations, there are various technical issues underlying the process. One of these resides 
in making accurate measurements through credibly small error bars. Moreover, one needs to 
have contact with the continuum structure of the same quantities once the lattice regularization 
is lifted. Indeed it should be the case that when low energy estimates from the lattice are 
extrapolated to the high energy continuum evaluation there should be reasonable agreement. 

To this end over a period of years there has been ongoing comparison of lattice computations 
with high energy perturbative expressions. The latter are determined in the chiral limit, usually 
in the MS scheme, to as high a loop order as is calculationally possible. Indeed with the current 
loop calculation technology this invariably means three loops through the use of the Mincer 
package, [3 [8]. As highlights of this bridge we mention quark current renormalization for zero 
momentum operator insertions, p2 [10], and more recently at non-zero momentum insertion, 
[1 H 112 1 [T3 l 114 1 115]. The latter are known at two loops at the symmetric subtraction point for the 
MS scheme. However, in both cases results have been determined in the regularization invariant 
(RI), PHUT], class of schemes which includes RF and RF/SMOM, [TT], where SMOM indicates 
momentum subtraction at the symmetric point. These are schemes which are devised for lattice 
regularization in order to minimize the use of derivatives, which is computationally intense, in 
extracting renormalization constants as well as amplitudes. Both sets of schemes have continuum 
analogues, [HI IIP] . Building on this project it is natural now to turn from quark-antiquark 
operators representing mesonic physics to 3-quark operators for baryon problems. Therefore, it 
is the purpose of this article to compute the Green's function for the zero momentum insertion of 
the 3-quark operator representing a proton in the chiral limit to two loops in the MS scheme. In 
addition to assist with running to high energy we will compute the operator anomalous dimension 
to three loops in the same scheme. Both results will therefore be important for matching lattice 
results of the same Green's function in the high energy limit. For the operators which we 
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consider here, we restrict ourselves to the MS scheme and will not introduce any variant of the 
regularization invariant class of renormalization schemes. Though from the results compiled 
here it is possible to define such schemes for the 3-quark operator renormalization and then 
convert to the MS scheme. However, one reason for not choosing to explicitly include, say, RI' 
results rests in the fact that as noted in [14] the definition of such schemes is not unique for 
Lorentz structure beyond the simplest quark mass operator. Indeed in the tensor current case 
a scheme alternative, [14], to the original one of [11} [I~3] appeared to converge quicker at two 
loops. Only a three loop calculation would give more insight into this. 

In more detail we will evaluate the two loop Green's function for an arbitrary linear covariant 
gauge at a point where the squared momenta of the three external quark legs are all equal 
to the same non-zero value. Hence, there is no zero momentum external quark leg. This is 
appropriate since a zero momentum quark would be difficult to incorporate on the lattice. Several 
lattice studies of 3-quark operators and their low moments for similar Green's functions deserve 
mention. In [18] an initial one loop lattice analysis was performed motivated by proton decay in 
an SU(5) grand unified theory. More recently the QCDSF collaboration extended aspects of that 
analysis to a full examination of the 3-quark lattice operators including moments, [19]. Though 
in both cases precise measurements require continuum perturbation theory for matching. While 
the lattice ultimately will only require results in the Landau gauge, we choose an arbitrary 
gauge for internal checking purposes. For instance, in extracting the anomalous dimension of 
the inserted 3-quark operator, the result has to be independent of the gauge parameter in the 
MS scheme. This will represent a useful check here. Previously the two loop MS anomalous 
dimensions were computed in the Feynman gauge in [20]. The initial one loop analysis was 
carried out in [lj El [3] . In finding total agreement with the expressions of [U El H9 [20] at two 
loops we will have another check on our results. At three loops to reduce the computation 
time for the large number of Feynman graphs to be evaluated, we will restrict the calculation 
at that order to the Feynman gauge. Though the way that part proceeds there will be an 
internal check over and above that provided by the renormalization group equation. It is also 
worth noting that the anomalous dimensions of 3-quark operators have been determined to two 
loops in heavy quark effective theory, |21[ [22] . In that analysis one of the three quarks in the 
operator is regarded as having a mass significantly larger than the other two. Equally 3-quark 
operators have been used to estimate baryon masses using the operator product expansion, 

[231 EH ESI ESI EH EH ES] . 

While this summary covers our aims it would be remiss at this stage not to mention several 
technical problems which have to be addressed. The 3-quark operators share a similar feature 
to four-fermi operators. Not only can they mix under renormalization but within dimensional 
regularization in d = 4 — 2e spacetime dimensions, which we use throughout and where e is 
the regularizing parameter, the operators will mix into evanescent operators. These exist in the 
analytically continued spacetime but are non-existent in strictly four dimensions. However, their 
presence within the renormalization has to be taken into account. We will use the projection 
technique of [30]. Coupled to this is the underlying 7 s issue due to the fermion handedness. For 
automatic symbolic manipulation calculations a method was developed in [31J to accommodate 
7 5 with operator renormalization. We adopt and adapt that technique for the 3-quark operators 
relating to the proton. At two loops the 7 s problem does not arise in MS for the operator 
anomalous dimension due to the nature of the mixing matrix. At three loops the problem will 
be evident and needs to be treated. 

Finally, we will extend a more recent two loop renormalization of a 3-quark operator which 
was introduced in [32]. In [32] rather than use an initial operator with the correct quantum 
numbers the most general 3-quark operator, devoid of any 7-matrix structure, was renormalized 
using dimensional regularization and an MS scheme subtraction. Clearly the seed operator 
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will mix into a set of 3-quark operators with 7-matrices. These arise from the internal quark 
propagators and quark-gluon vertices of the Feynman diagrams at each order. The claim in [32] 
is that from this mixing one can overcome not only 7 s issues but also handle the evanescent 
operator problem in a systematic way. This is motivated by earlier work of [33]. Hence we will 
first check the two loop general anomalous dimension of [32] in an arbitrary linear covariant 
gauge and then extend the result of [32] to three loops. The former will also play the role of a 
subsidiary check on our proton operator MS result. Though it is worth stressing that the scheme 
dependent parts of the two loop results of [32J are not the same as those of [20| which is regarded 
as being the MS scheme. As noted in [32] they have provided a conversion factor which appears 
to derive from O(e) corrections similar to contributions from evanescent terms. However, one 
advantage of the general operator of [32J is that the anomalous dimension of 3-quark operators 
of the spin-(j,j) Lorentz type can be easily deduced from the general anomalous dimension. 
Hence we will provide the three loop value for, say, the Ioffe current [23] among other quantities. 

The article is organized as follows. Section 2 is devoted to the background for the renormal- 
ization of the operators we will examine at three loops including the definition of the operator 
basis in (^-dimensions. We provide technical details of the calculation of the Green's functions 
of interest in both momentum configurations in section 3. The results of the renormalization 
are recorded in section 4 including the finite renormalization required to treat 7 s . The explicit 
form of the two loop amplitude relevant for lattice matching is presented in section 5. Section 
6 is devoted to the extension of [32] to three loops where we record the anomalous dimensions 
of various spin operators at three loops in the scheme used in [32]. Finally, we conclude with 
remarks in section 7. 



2 Background. 

We begin by first discussing the actual operators we will renormalize. For further background 
we note that we have based this on the analysis of [20J . Clearly the proton involves up and down 
quarks and the relevant operators in the chiral limit with the correct parity are 

udu = e IJK u I ftjyT Cj 5 d Kj (21 ) 



which carry one free spinor index. Here C is the charge conjugation matrix which satisfies 
CC = - 1 with 

C(^fC = 7^ . (2.2) 

The indices /, J and K are SU(3) colour indices and throughout we work in this specific 
Lie group. Therefore, within our calculations products of the group generators are simplified 
automatically with 
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(2.3) 



where 1 < a < 8. This means that when, for example, we make use of the quark wave function 
and gauge parameter anomalous dimensions as well as the QCD /3-function then the usual group 
Casimirs, Cf, Ca and Tp, of those expressions are also evaluated at their SU(3) values. Indeed 
we stress that the 3-quark operators we consider here have no physical meaning for colour groups 
other than SU(3) since only for that group are they gauge invariant. 



In (|2.ip we have two operators of the same dimension and discrete symmetries. These will mix 
under renormalization. However, various linear combinations will produce the correct number of 
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left and right handed quarks to produce the operator which actually corresponds to the proton 
itself. As was discussed in [20] for practical computational purposes it is more appropriate to 
use a related set of operators to perform our two and three loop calculations. These operators 
are given by 

O x = e /Ji V (Vf 'C^ K ) 

2 = e IJK ^ (i^ J ) T C^ K ) , (2.4) 

where we have omitted the flavour indices, and are related to the operators of (|2.ip by 

d = 7 5 Of u . (2.5) 

From now on our focus will be on Oj. Though we note at this point that whichever set of 
operators one uses they are gauge invariant for SU (3) colour and so in the MS scheme, which 
we will use, the anomalous dimensions will be independent of the gauge parameter. This will 
be checked explicitly to two loops. The original computations of [20] were performed in the 
Feynman gauge rather than the arbitrary linear covariant gauge we use here to two loops. 

For instance, as we have two operators of the same dimension and symmetries their mixing 
under renormalization is handled by a mixing matrix of renormalization constants and thence a 
matrix of anomalous dimensions. For the former we therefore have 



O lo = Z l3 O j (2.6) 

where the subscript denotes the bare operator. In our conventions the matrix of anomalous 
dimensions is defined by 

lijia) = - fi—lnZij (2.7) 

where 

d d d 

fi— = (3(a)— + a-y a (a,a)— (2.8) 
d\i oa oa 

and li is the renormalization scale introduced to ensure that the coupling constant, g, is dimen- 
sionless in d-dimensions. Our gauge parameter is a with a = corresponding to the Landau 
gauge. We choose to work with the coupling constant a which is related to the gauge coupling 
constant and the strong force coupling constant, a s , by 

(j 2 a, . 

(2.9) 



16vr 2 ' 4vr ' 

For mass dependent renormalization schemes 7y(a) can depend on a but we have omitted any 
a dependence on the left side of (|2,7p since we focus on MS. For practical purposes in the 
extraction of the anomalous dimensions we note that (|2.7|) implies 

1ik(a)Z k j = Z ik ^ kj (a) = - ^-^ z ij (2.10) 

and we have checked that the same mixing matrix emerges irrespective of which way the matrices 
are multiplied. 

It should be stressed that this discussion about the structure of the mixing matrix, where i 
and j run over 1 and 2, is in essence the situation in four spacetime dimensions. However, as 
will be apparent later since we will be using dimensional regularization the mixing matrix is not 
finite dimensional in (i-dimensions. It will have to be extended to an infinite dimensional case 
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at all orders in perturbation theory. This is because in d-dimensions operators will be generated 
through the renormalization which have no strictly four dimensional equivalent. Put another 
way their four dimensional equivalent is the zero operator and so such operators in d-dimensions 
are referred to as being evanescent. Their presence in this analysis is not solely because of the 
presence of 7 5 . These evanescent operators would arise if there was no mixing and one was 
only considering 0\. However, both the evanescent and 7 s issues will have to be dealt with 
at the appropriate point. While we have noted that Z^j will become infinite dimensional in 
d-dimensions, it will do so in a controlled way in that at each order in perturbation theory the 
matrix will be extended by the appearance of a new operator. In other words at each order in 
perturbation theory the matrix enlarges but remains finite dimensional. 



To appreciate these remarks we need to focus on the treatment of 7-matrices. In d-dimensions 
one has to extend the basis of 7-matrices to an infinite set of matrices denoted by rj^C" , 
|30[ I3"4"l [35] , which are totally antisymmetric in the Lorentz indices and defined by 

T ni n n = 7 Un (2.11) 

where a factor of 1/n! is understood and n is an integer, n > 0. These generalized matrices 
span spinor space in d-dimensions and the underlying algebra necessary for loop calculations 
has been developed in various articles, [301 EH ES] . For instance, the trace operation is isotropic 
with respect the basis since, [361 [37], 

tr fev^r^-M oc 5 mn I^-^-^ . (2.12) 



(m) (n) 

\T rnn /" \ v l (Vi Tif\l tyi nrvi/inp f~i c r\ ~r~i ~vt it/i cn v v ^ /-yt rnv I 

(n) 



Moreover, it is possible to write products of the original 7-matrices as a finite sum over r^ 1 " / '" 



This can be achieved recursively by applying the relations, [351 1361 13"T] . 



n 



(n) ' (n+l) / .A i ' (n— 1) v ' ' 



r=l 
n 



' (n) (n+1) / / ' (n— 1) » ' / 



r=l 



where 77^^ is the spacetime metric tensor. For reference when one restricts the -matrices to 
four dimensions we have 



L (2) 



d=4 ~ * ' W 



d=4 



= for n > 5 (2.15) 

d=4 



where e Mycrp is the four dimensional totally antisymmetric pseudotensor. As a note we mention 
that the 7 s matrix which exists in strictly four dimensions and defines chirality has absolutely 
no connection whatsoever with p^—w m d-dimensions. The former object exists only in four 
dimensions and Larin's procedure, [31], which we use to handle 7 s within dimensional regular- 
ization will be detailed later for the present computation. As a notational comment we will use 
7^ and r^s synonymously in d-dimensions since in (|2.14j) as the former is less clumsy but regard 
as the purely four dimensional object. 

There are two sources of these generalized 7-matrices within the computations. The first is 
the simplest and alluded to already and that is that the product of 7-matrices which remain in 
the calculation of either Green's function can be written in the general basis. This leads to the 
second source which is the generation of operators into which the seed operator, (|2.4p . of the 
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Green's function mixes under renormalization. In d-dimensions given that we have to use the 
generalized basis we then have to use the generalization of the four dimensional operators to the 
same basis. Therefore we define the new (i-dimensional operators 

0(n) = ^(r5"""/)((^) T Cr (n)w ... /In /) , forn/4 

0(4) = ^^(rf 4 f^ 4 ^)(M T Cr (4)w ^ 4 ^) , forn = 4 (2.16) 

where the factor for is to ensure that there is a direct mapping in the four dimensional limit 
to the original operators. This follows from (|2.15p and the fact that in e^ uup e^ uap = 24. Hence, 

( 4) = 2 + O(e) (2.17) 

in the limit to four dimensions. As a mnemonic we note that 

O n = C( 4n _4) (2-18) 



for n > 1 and 



02\ d=A = e 



O n \ d=A = for n > 3 . (2.19) 

Thus at the outset we are forced to consider a more general renormalization from which the 
anomalous dimensions of the original operators will emerge as a corollary. So within all our 
computations we will write the Green's functions in terms of T^-matrices. For the finite parts, 
which is ultimately required for lattice matching, the Lorentz indices on these matrices can 
be contracted with the external momenta. Though such objects will not be associated with a 
divergence in e since this would violate the renormalizability of the operators. Before restricting 
the (i-dimensional Green's function there are more channels. This is one of the reasons why 
we have not taken the projection approach which was the main tool for mesonic operators, 
115] . For instance, one has to have knowledge of the full basis to say three loops for (|2.26p 
and then construct the projection tensor to isolate all the possible spinor channels. Within a 
symbolic manipulation approach this would significantly reduce run times due, in part, to having 
to internally manipulate products of r( n )-matrices. Instead we have constructed the relations 
between the products of rj^y Mm and T^ Mn for various values of m and n which occur. These 
have been encoded within a module in the symbolic manipulation language, Form, [38J which we 
use throughout. Such a product can be written in terms of r^"' <Jp where \m — n\ < p < (m + n) 

and Cj G {fJ-i, . . . , fjL m , v\, . . . , u n }. The indices not used from this set in " <Jp appear in the 
rj^ u tensors which are required to keep the total number of free indices of each term as (m + n). 

Given the appearance now of the full basis of operators Or n \ in ci-dimensions we will de- 
termine the associated renormalization constant matrix Z^. This matrix will increase in size 
at each loop order but from it we will determine what we will refer to as the naive anomalous 
dimension matrix, jij (a) . It would ordinarily correspond to the correct four dimensional anoma- 
lous dimensions but in using these generalized operators we have ignored the problem of 7 5 as 
well as the effect of the evanescent operators. The latter affect the structure of the four di- 
mensional anomalous dimensions even though the evanescent operators are non-existent in four 
spacetime dimensions. To account for this we note the formalism developed in [30] which ap- 
pends to the naive anomalous dimensions extra contributions which derive from the evanescent 
parts. We make minimal comment on the technique here since it transpires that the effect they 
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have on the proton operator renormalization will not occur until four loops. Though if there 
were contributions these together with the naive anomalous dimensions would contribute to the 
correct four dimensional result. We say contribute as one has also to deal with the absence of 
chiral symmetry in d-dimensions which has been ignored with the choice of . Including the 
evanescent effects with the naive anomalous dimensions for the present case would not produce 
a result consistent with chiral symmetry in strictly four dimensions. The procedure we have 
chosen to do this is based on Larin's method, |31j . which was developed for flavour non-singlet 
and singlet quark currents as well as the chiral anomaly. We will discuss the technical aspects 
of the calculation for our case later but in essence one needs to append a finite renormalization 
constant to the naive renormalization constant. It is chosen in such a way that in strictly four 
dimensions the anti-commutativity of 7 s with 7^ is restored. This aspect is treated after the 
contributions from the evanescent operators have been included. The criterion for defining the 
condition will be similar for our operators Oi, i = 1 and 2, and like |31j will be derived from 
knowledge of the finite part of the Green's functions with the operators inserted. Though the 
finite renormalization can also be derived from the difference in anomalous dimensions. How- 
ever, in either case unlike [31] it will be a matrix of finite renormalization constants leading 
to an additional matrix of anomalous dimensions, 75,ij( c 0- Thus the correct four dimensional 
anomalous dimension matrix for (12.11) or (12. 4p will formally be 

Hj{a) = %(a) + 75,i.7 (a) (2.20) 
where we exclude any evanescent part for this case. Here 

I5,ij(a) = - li^-Z% (2.21) 

and Zfj is the finite renormalization constant matrix. If the original seed operators had included 
additional 7-matrices, with or without free Lorentz indices such as those of (24, 39], then the 
evanescent operator contribution could occur at the three loop order we are interested in here. 
In writing 75^ (a) we are assuming that the finite renormalization is independent of the gauge 
parameter. In the cases examined by Larin in |31j and in a more recent analysis of diquark 
operators, [3D], the finite renormalization did not depend on a which is what we found here. 
Indeed it is worth noting that as this additional piece corresponds to a finite renormalization, 
for the three loop mixing matrix in four dimensions one only requires Zfj to two loops because 
of the presence of /3(a) in (|2.2ip . 

We close this section by summarizing the two calculations we carry out. First, in each case 
a 3-quark operator is inserted at zero momentum into a quark 3-point function. For illustration 
if for the moment we denote this generic operator by O then we will calculate 

(Mp)MQ)^i r )°s(0)) (2.22) 

to two loops to the finite part for arbitrary gauge parameter a in the MS scheme. This is to 
assist with lattice matching to the same quantity in the Landau gauge. For this case there are 
3 one loop and 40 two loop Feynman graphs to compute and we will call this the symmetric 
setup. In (|2.22p momentum conservation implies 

r = - p - q (2.23) 

and we use a symmetric subtraction point for the external legs 

p 2 = q 2 = r 2 = - fj, 2 (2.24) 
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which implies 



1 2 

PQ = 2^ • 



(2.25) 



We have chosen the subtraction point to be (—/J 2 ) so as to omit logarithms in the finite part of 
the Green's function. In order to extract the three loop mixing matrix of anomalous dimensions 
we have to consider a different momentum configuration so that the Mincer algorithm, [7J, can 
be applied. We will refer to this as the Mincer setup. In this case the Green's function is 



and there are 784 three loop diagrams to determine in addition to those noted earlier at lower 
order. In both situations, (I2.22P and (|2.26p . the Feynman diagrams are generated automatically 
using the Qgraf package, jH], before being converted into Form, [55], input notation. The 
latter is the symbolic manipulation language we use and the full automatic computation is 
written in terms of it. Though the compilation of the expressions for each Feynman diagram for 
both setups were run with the threaded version, Tform, |42| . The nullification of an external 
quark leg momentum in effect produces a 2-point function which Mincer requires. 

3 Computational technicalities. 

We devote this section to the various technical aspects of the computation. In order to evaluate 
the Green's function for the full momentum or symmetric configuration we use the Laporta 
algorithm, |43| . Each of the one and two loop Feynman graphs we have to evaluate involve 
strings of 7-matrices which have external and internal momenta embedded within them. We 
have proceeded by writing all integrals in terms of scalar integrals. By this we mean integrals 
where there are at most scalar products of the momenta and the strings of 7-matrices have 
contractions with only external momenta. In other words these only play a passive role in the 
subsequent evaluation and are written in terms of the generalized r( n )-matrices. To achieve this 
we make a projection of the tensor integrals onto a basis of tensors built from 77^, and q^. 
For the tensor reduction we use we have at most a rank 5 tensor built from the two internal 
loop momenta at two loops. This is due to the fact that we are computing for an arbitrary 
linear covariant gauge. A Feynman gauge computation would be more compact but would have 
limited applicability for the lattice. Once the scalarized integrals have been determined these 
are rewritten purely in terms of the propagators of the graph and any additional propagators 
which are not part of the topology, [43] . These latter propagators are required for irreducible 
scalar products but are chosen in such a way as to cover all possible scalar products of the 
internal momenta with themselves and the external momenta. In this form one applies the 
Laporta reduction. This is an algorithm which systematically constructs all the integration by 
parts relations and optionally the Lorentz identities between all the integrals which are needed. 
From this tower of relations it is possible to algebraically relate all the scalar integrals to a base 
set of master integrals. These are evaluated by explicit integration and thus the evaluation of 
the Feynman graph is complete. 

In describing the general procedure we note that for practical purposes one has to use 
computers to implement the Laporta algorithm. We have used Reduze, [44], which uses the 
GiNaC computer algebra system, |45] . and is written in C++. For the one and two loop graphs 
we need to evaluate it transpires that there are three basic topologies. There is one at one loop 
and two at two loops. For the latter one is the ladder graph and the other is the non-planar two 
loop 3-point graph. These and their extension to include one other propagator are sufficient to 
cover all possible irreducible tensor integrals. Thus using the Reduze package we have created 



p2 — _p2 



(2.26) 
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a database of relations covering all possible levels of scalar integrals which can arise. The ones 
which are explicitly required are extracted and converted into a Form module which is called at 
the appropriate point of the automatic computation. At the end the explicit expressions for the 
master integrals are substituted. The ones we use are summarized in [13] but were evaluated in 
various articles, [46j HTJ HSJ [49], using a variety of techniques. 

For the Mincer situation, (|2.26p . the momentum configuration with one nullified external 
quark leg could potentially introduce spurious infrared infinities which would need infrared 
rearrangement. However, this does not arise. This would be the case if the Feynman integral 
produced propagators such as l/(fc 2 ) 2 , where k is an internal loop momentum, which are infrared 
singular. These are absent because the quark propagator retains ^ in the numerator or the triple 
gluon or ghost vertices carry a momentum to lift the potential infrared singularity. By contrast 
to (|2.22j) we take a more general operator in order to compute the anomalous dimension. In 
particular we seed the Green's function with 



which is not decorated with 7- matrices and a, j3 and 7 are spinor indices. There are several 
reasons for doing this. One is a practical one to do with the size of the three loop calculation. 
There are 784 three loop diagrams to determine and thus to keep computer run times to a 
minimum it transpires that it is more efficient to evaluate the diagrams with (|3.ip and then 
introduce the 7-matrix structure appropriate to each original operator of (|2.16p when summing 
the diagrams. In addition it also allows one to quickly construct that part of the mixing matrix 
relating to the evanescent operators which are generated rather than have to repeat a full run. 
Indeed in that case the calculation would necessarily be slower as it could involve rj^j"^' 1 for 
values of n up to 20 for each of the three loop graphs. Another reason for proceeding with (13. II) 
is that we can extend the recent calculation of [32] at the same time. 

However, in choosing to calculate in this more general way there are several technical issues 
to be overcome which are rooted in the Mincer algorithm. It computes massless scalar 2-point 
functions to three loops in dimensional regularization. With the 3-quark operator zero momen- 
tum insertion in a quark 3-point function the nullification of an external quark leg momentum 
ensures we have in effect a 2-point function immediately. To obtain scalar integrals within the 
symbolic representation of each Feynman diagram we strip off the 7-matrix structure from the 
numerator as for (|2.22p . An alternative approach would be to project out the Lorentz structure 
but this is too cumbersome especially as one has to have prior knowledge of the full structure of 
the Green's function at each loop order. Removing the 7-matrices instead leaves each diagram 
as a sum of Lorentz tensor integrals. Given that we are in the chiral limit these integrals will 
be of even rank. In the Feynman gauge they will be rank (21) where I is the loop order. By 
contrast for a general linear covariant gauge these integrals will be at most rank (4/). In addition 
the Mincer algorithm, [7], codes the internal loop momenta for each of the Mincer topologies 
in terms of its own labelling. In particular each line of a topology is asigned an internal mo- 
mentum label pi. The conservation of energy momentum at each vertex is then encoded within 
the integration routine for that topology. In other words there are no Mincer labels such as 
(Pi ~ P2) or (Q — P3 — Pe) where Q is the external momentum of the 2-point function. Thus all 
the tensor integrals involve products of internal momentum vectors pi where 1 < i < 8 for three 
loops. For lower loop orders there are fewer internal momentum labels. Therefore, the problem 
of evaluating each graph of the Green's function requires converting these even rank tensors into 
scalar integrals, which are straightforward to compute in Mincer, and Lorentz tensors built 
from r)uv an d Q^- The procedure for this is straightforward. Using only knowledge of the rank 
we have written down the most general tensor basis for each rank and then determined the 
scalar integral amplitude by the method of projection for each tensor. For instance, for rank 



O = e IJK ^U, 
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2, 4 and 6 there are respectively 2, 10 and 76 tensors. Using Form we have constructed the 
decomposition for an arbitrary numerator which is straightforward. For rank 8 the projection 
matrix would be 764 x 764. Rather than use a decomposition into r)n V and we chose the 
transverse and longitudinal projection tensors 

PAQ) = *ir-^r > VW) = ^ (3.2) 

which satisfy the simple properties 

PAQ) + L AQ) = v»v , PAQ) pua (Q) = p „ a (Q) 

P^(Q)L™(Q) = , L^(Q)L™(Q) = L/(Q) . (3.3) 

The benefit of this is that the full matrix becomes block diagonal with the largest submatrix be- 
ing 105 x 105. Once the general decomposition has been derived it is encoded as an integration 
module in FORM at the appropriate point in the overall algorithm. With rank 8 decomposi- 
tion the anomalous dimension matrix can be determined to two loops in an arbitrary linear 
covariant gauge. At this point we note that since the operators we consider are gauge invariant 
this allows us to check that in MS the two loop anomalous dimensions are independent of the 
gauge parameter. This is a strong check on the derivation of the Form module encoding the 
mapping of the Lorentz tensor integrals to the scalar integrals. Indeed the two loop result of 
|20j was performed only in the Feynman gauge. For three loops one could in principle extend 
the construction to rank 12. Instead we have taken a different approach and chosen to calcu- 
late the three loop diagrams solely in the Feynman gauge. In this case this only requires the 
decomposition up to rank 6 which has already been tested for the arbitrary gauge calculation 
at two loops. Moreover, in the Feynman gauge the diagrams can be evaluated more quickly. 
The independent check on the three loop anomalous dimensions will be consistency with the 
renormalization group equation as the double and triple poles in e are already fixed from the 
one and two loop renormalization constants. We believe this calculational approach from the 
point of view of exploiting the properties of the gluon propagator in various gauges is the most 
efficient to deduce the anomalous dimensions. However, the finite part of this Green's function 
at three loops is not useful for lattice matching as the Landau gauge results will not be deter- 
mined. Indeed it is not clear whether they would be meaningful anyway in this instance as the 
nullification of the momentum of an external quark leg could be problematic from the point of 
view of infrared divergences in the chiral limit. 



4 Three loop anomalous dimensions. 

While we have split the discussion on the technical aspects of each setup we collect the results 
of our renormalization in one section. This is because for both calculations we obtain the 
same results to two loops for an arbitrary linear covariant gauge. As indicated in the Mincer 
discussion we have reasonable internal checks on the three loop part of the anomalous dimensions 
from the renormalization group equations. First, we record the naive anomalous dimensions for 
the seed operators Or a \ and Ou\ where the labels 1 and 2 refer to these operators respectively. 
To three loops in MS we have 

a 2 

7U (a) = -2a-[2N f + 51] — 

y 

a 3 

+ [260AT/ + [4320C(3) - 4656]^/ + 1296C(3) + 23481] + 0(a 4 ) 
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1 fl n 3 

712(a) = - a 2 + [216C(3) - 153 - UN f ] - + 0(a 4 ) 

1 D n 3 

721(a) = ya 2 + 4 [18C(3) + 331 - 22N f ] - + 0(a 4 ) 

a 2 

7 22 (a) = -2o-[2JV / + 51] — 



a 3 



+ [260AT/ + [4320C(3) + 1344]iV) + 1296C(3) - 75519] — + 0(a 4 ) (4.1) 

where C,{z) is the Riemann ^-function. The two loop parts are calculated for arbitrary a and 
agree with the full two loop mixing matrix of [20J. This is because at this order the 7 s problem 
can be ignored as was noted in [2D] and hence the two loop naive anomalous dimensions are 
sufficient to determine the proton operator wave function renormalization. As we generate 0( 8 ) 
or O3 we have to include the next row and column of Jij(a) to the appropriate orders in a. Thus 
we have 

a 2 

733(a) = - 2a - [2N f - 19869] — + 0(a 3 ) 

723(a) = ^a 2 + [4824C(3) - 2745 - 147V;] ^ + 0(a 4 ) 
73i(a) = 732(a) = 0(a 3 ) 

713(a) = 0(a 4 ) . (4.2) 

The key terms in this are those in the upper triangle. In the formalism of [30J 731(a) and 732(a) 
correspond to the generalized /3-functions used in the evanescent projection formalism to extract 
the effect the evanescent operators have on the strictly four dimensional anomalous dimensions, 
before we restore chiral symmetry and the anti-commutativity of 7 s . As these two anomalous 
dimensions are zero to 0(a 3 ) then there can be no evanescent operator affect until four loops in 
the full four dimensional anomalous dimensions. 

In order to complete the computation we require the finite renormalization constants which 
restore the anti-commutativity of 7 s in strictly four dimensions. As background it is worthwhile 
recalling the Larin procedure, [31], for extracting the finite renormalization constant for the 
flavour non-singlet pseudo-scalar quark current, O m5 = i/j^i/j. Here we assume the operator is 
inserted in a quark 2-point function at zero momentum. With a naive anti-commuting 7 s in 
c?-dimensions one would simply anti-commute the 7 s out of all the diagrams. However, retaining 
it inside diagrams and using the generalized FV n ) -matrices to compute the Green's function 
one obtains a finite Green's function in the four dimensional limit after the renormalization 
constants have been determined in the MS scheme. The value for this should be equivalent to 
that calculated with the naive anti-commuting 7 s in d-dimensions. However, it is not the same 
and therefore to remove the discrepancy one defines a finite renormalization constant, |31] . If 
we define this as Z 5 for the pseudo-mass operator then more specifically the condition is defined 
by 

7 5 ^(p)o m m(-p))\ d=4 = z 5 (ip( P )o m5 m(-p))\ d=4 (4.3) 

where O m = ipip is the quark mass operator. While this is the situation for a single operator 
with no mixing the renormalization of the 3-quark operator is complicated by the mixing. In 
this respect it is notationally difficult to extend the definition of (|4.3p to the matrix case. At a 
formal level the definition can be represented by 

(Oi) = Z 5 ( 7 5 ® 7 5 ) (O2) (4.4) 

where the tensor product of the 7 s matrices acts on the appropriate spinor indices of the Green's 
function. Also Z 5 represents the finite renormalization matrix and its effect within the finite 
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part of the Green's function which is the meaning of the multiplication. While this is the 
normal Larin procedure where the naive MS anomalous dimensions are computed first and then 
the finite renormalization, in practical terms it is quicker to do both processes together. In other 
words we absorb a certain finite part into our MS renormalization constants. The condition for 
this is that the finite parts of both Green's functions after renormalization in four dimensions 
are equivalent up to multiplying by 7 s (g> 7 s . However, to assist others who wish to reproduce 
the results we will present them as in the two stage Larin method. Thus in addition to the naive 
anomalous dimensions we have extracted the anomalous dimensions associated with the finite 
7 5 renormalization matrix. These are 



where the linear factor in Nf derives from the one loop /3- function for SU(3). In defining the 
finite renormalization condition (|4.4p we have not included the momentum configuration of the 
external quark legs. This is because (j4.5j) has been derived for both setups. The emergence of 
the same finite renormalization from both calculations is an independent check. Moreover both 
computations were performed in an arbitrary linear covariant gauge. However, as is evident 
from fj4.5j) the result is independent of a. For all the quark bilinear current operators considered 
in [31] and the diquark operators examined in |40j which contained 7 s , the associated finite 
renormalization constant was also independent of the gauge parameter. Indeed this follows from 
a simple observation that the finite renormalization constant in the bilinear current cases is 
related via renormalization group arguments to the naive anomalous dimensions of the particular 
operators involved in the defining condition. In the case of our earlier example this would be 
the naive anomalous dimensions of O m and O m5 . As the operators are gauge invariant and 
being renormalized in a mass independent scheme then as the operator anomalous dimensions 
are independent of the gauge parameter then so to is the finite renormalization constant. By 
contrast if the original operator renormalization had been performed in a mass dependent scheme 
then not only would the operator anomalous dimensions be gauge dependent but the finite 
renormalization constant would too. This is of course barring accidental cancellations of the 
gauge parameter which is in principle possible. For the case of the operators 0\ and O2 a 
similar argument can be established. Though due to the mixing it derives from within the formal 
renormalization condition (I4.4p and the specific way we have defined the finite renormalization 
constant within the explicit calculation of the Green's function of each operator in the one stage 
method of applying the Larin technique which we used here. To summarize the upshot of that 
analysis translates into the equations 



These equations are analogous to the situation in the quark current example. As each of the 
naive operator anomalous dimensions are gauge independent, since the operators themselves 
are gauge invariant, and we are in the MS scheme, which is a mass independent scheme, then 
the two finite renormalization constants are automatically also gauge independent. Indeed from 
the explicit expressions for the naive three loop anomalous dimensions it is straightforward 
to check that (|4.5p are consistent with these general expressions. This is a reassuring check 
since we derived (I4.5P from the finite part of (|4.4p which was calculated in an arbitrary linear 



75,21 (a) 



75,11 (a) 



75,i2(a) = 0(a 4 ) 

125 [2^-33]^ + 0(a 4 ) 




(4.5) 




(4.6) 
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covariant gauge in both calculational setups to two loops. It was only the three loop Mincer 
calculational which was carried out in the Feynman gauge. Therefore, there is consistency with 
that calculation too. 

Equipped with this finite renormalization we can now determine the full four dimensional 
mixing matrix using the naive anomalous dimensions. We find to three loops that 

. r 2 

- 2a- 



7n(a) = 722(a) 



[2AT/ + 51]- 

+ [260-/V/ + [4320C(3) - 4656]^/ + 1296C(3) + 23481] 
+ 0(a 4 ) 



a 

162 



712(a) = 721(a) 



10 



a 2 + [216C(3) - 153 



UN f ] - + 0(a 4 ) 



(4.7) 



The effect of the finite renormalization has been to restore the symmetry of the mixing matrix 

so that the diagonal entries are equal and the off-diagonal are the same but different to the other 

two. This structure was present at one and two loops, [U [21 [31 [20]. However, for the actual 

proton anomalous dimension we need to have the correct handedness of the up and down quarks 

which requires the eigen- anomalous dimensions which are 

2 3 

7+ (a) = - 2a - [2N f + 21] — + [260iV f 2 + [4320C(3) - 4740]iV f + 2592C(3) + 22563] — 

9 J 162 

+ 0(a 4 ) 



7- (a) 



where 



3 

- 2a - [2N f + 81] ^- + [260iV/ + [4320C(3) - 4572]A> + 24399] ^— 
+ 0(a 4 ) 



(4.8) 
(4.9) 



7+ (a) = 711(a) + 712(a) , 7_ (a) = 711(a) - 712(a) . 

It is the latter, 7_(a), which corresponds to the proton. Numerically we have 

7+ (o) = - 2.0000000a - [0.2222222iV/ + 2.3333333] a 2 

+ [1.60493837V/ + 2.7955915iV) + 158.5106882] a 3 + 0(a 4 ) 

7_(o) = - 2.0000000a - [0.2222222iV/ + 9.0000000] a 2 

+ [1.60493837V/ + 3.8326285^ + 150.6111111] a 3 + 0(a 4 ) . 

Having established the anomalous dimensions we can construct the renormalization group in- 
variant current, j_, using the same notation as |20| . It is defined in the conventional way by 



(4.10) 



cxp 



7-(z) 



and satisfies 



Solving 



dj_ 

dfj, 



7i-(o)i-(M) 



(4.H) 

(4.12) 
(4-13) 



explicitly we formally have 



Ij- (a) 



,71 //3i 



1 + 



11 



Pi 

73 _ 
Pi 01 
+ 0(a 3 )] 



+ 



71A t| 
+ /?i 2 + 



PI 



71^ 
Pi 



27i72/?2 + llPl 



Pi 



PI 



a 

T 



(4-14) 
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where 



7_(o) 
f3(a) 



7i a + 72a 2 + 73« 3 + 0{a A ) 
(3 ia 2 + (3 2 a 3 + /3 3 a 4 + <3(a 5 ) 



(4.15) 



This produces 
7j- («) = 



1 -[4^-588^+ 2835] 3|2iV/ _ 33p 

+ [2080AT/ + [34560C(3) - 157368]iV/ + [4596912 - 1710720C(3)]JV/ 
+ [28226880C(3) - 70113330]iV/ + [580876920 - 155247840C(3)]iV / 



- 1825381251] 
For three flavours this gives 



+ 0(a 3 



77- 



\N f =3 



1 



41 
— ( 
81 



108[2iV) - 33] 4 

[116640C(3) + 275215] 2 



6/I33-2JV/ 



26244 



a 2 + 0{a A 



,2/9 



or 



75-(a)U r 



[1 - 0.5061728a - 15.8292531a 2 + 0(a 3 )] 



,2/9 



(4.16) 



(4.17) 



(4.18) 



numerically. To gauge the effects of the two loop correction we can compare the numerical value 



of 7j_(a) 



7V /= 3 



at one loop with that at two loops for a s = 0.1. We find that the two loop 



correction modifies the one loop value by around 0.1%. 



5 Amplitudes 



In this section we record the explicit values of the Green's function (|2.22|) to two loops in an 
arbitrary linear covariant gauge in the MS scheme. This represents one of the main results of 
the article as it will be of use for lattice matching. We havqj 



(^a{p)Ml)^-P-l)°ls(0)) 



symm 



1 + - [2vr 2 a + 15a + 2vr 2 + 15 



-3V'( 
+ [[2160^(i) 



- 3^'(i)a] a 

1440vr 2 - 19008] Nf + 1440^' (|) 2 



[8262a 2 + 23652a + 122742] ip' 



- 1920V>'(5)7r 2 + [711a + 207] ip" 
+ [136080a - 392688] s 2 (f ) 

+ [785376 - 272160a] s 2 (f ) 

+ [654480 - 226800a] s 3 (f ) 

+ [181440a - 523584] s 3 (|) + [88 - 1896a] vr 4 

+ [5508a 2 + 15768a + 81828] vr 2 

- [18468 + 12636a] S - [73224a + 23976] C(3) 



+ 188244a + 403461 + [2929 - 1015a] 



7f 

7! 



"The full analytic form of the amplitude for an arbitrary gauge and the anomalous dimensions, have been 
included in an attached electronic data file. 
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+ [32724 - 11340a] 



ln(3)vr 



+ [945a - 2727] 



ln 2 (3)7r 



2916 



+ 5 [[648C(3) - 1323 - 540V>'( ^vr 2 - 3888s 2 (§) 
+ 7776s 2 (§ ) + 6480s 3 (§) - 5184s 3 (§) 

+ 360vr 2 + 29^ + 324^ 
^3 y/3 

,ln 2 (3)7rl 



-27- 



a 
486 



+ 



4 [3ip'(l)a + 3^(|) - 2vr 2 a - 2vr 2 ] 
+ [[10567T 2 - 1584^(1)] N f 



81 



- [333a + 513] V (§) 

+ [2322a 2 - 1188a - 1350] ^'{\) 

- [104976a + 295488] s 2 (§) 
+ [209952a + 590976] s 2 (§) 
+ [174960a + 492480] s 3 (§) 

- [139968a + 393984] s 3 (§) 
+ [888a + 1368] vr 4 

+ [900 + 792a - 1548a 2 ] vr 2 

+ [4212a + 4212] E + [15552a - 5184] C(3) 



7T 



+ [2204 + 783a] -= 
V3 



+ [24624 + 8748a] 



ln(3)7r 

~7T 



[729a + 2052] 



ln 2 (3)7r 



7 2 1 l 



4374 



'a<5 



w ^ i > 

3' 



+ 5 [[648C(3) + 216V''( 1 r)vr 2 - 18*0 

- 3888s 2 (f ) + 7776s 2 (|) + 6480*-, 

-5184s 3 (f) + 48vr 4 

,ln(3)7r 



6 - 

144tt 2 + 29^ 



+ 324- 



27 



V3 
ln 2 (3)7r 



,2 "I 



4374 



7 a- 



/3 7 



+ 



2 [3^'(i)a + 3^(|) - 2vr 2 a - 2vr 2 ] — 

81 

+ [[912vr 2 - 1368^(1)] JV/ 
- [333a + 162] 
+ [2322a 2 - 1188a + 31050] 
+ [19440 - 104976a] s 2 (f ) 
+ [209952a - 38880] s 2 (f ) 
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+ 



+ [174960a - 32400] s 3 (§) 
+ [25920 - 139968a] s 3 (f ) 
+ [888a + 432] tt 4 + [4212a + 7128] S 
- [20700 - 792a + 1548a 2 ] vr 2 
+ [15552a - 19764] C(3) 
ln(3)7r 



+ [8748a - 1620] 



^/3 



+ [135 - 729a] 



ln 2 (3)7r 



+ [783a 



1451 W 



8748 
2vr 2 a - 



a 



p t 

ad 



a 



2vr 2 ]- 



81 



4[3V'(i)a + 3^(i) 

+ [[I968vr 2 - 2952^'(i)] N f 

- [666a + 765] <'(§) 

+ [4644a 2 - 2376a + 30780] 

- [295488 + 209952a] s 2 (§) 
+ [419904a + 590976] s 2 (f ) 
+ [349920a + 492480] s 3 (§) 

- [393984 + 279936a] s 3 (§) 

+ [1776a + 2040] vr 4 + [8424a + 11340] S 

- [20520 - 1584a + 3096a 2 ] vr 2 
+ [31104a - 21708] C(3) 

ln(3)7r 



+ [17496a + 24624] 



x/3 



[2052 + 1458a] 



ln 2 (3)7r 
V3 



IT" 



+ [2204 + 1566a] —= 
v3. 



8748 

+ [[[212^(1) - 144tt 2 ] N f + Ulip 



0", 



in f\\ 

3' 



31320^ (|) + 334368s 2 (f ) - 668736s 2 (f ) 
557280s 3 (f) +445824s 3 (|) - 1176vr 4 
20880tt 2 + 2916S - 17820C(3) 



- 2494 

+ 2322 
0(a 3 ) 



IT" 



V3 
ln 2 (3h 

~7T 



27864 



ln(3)7r 

~7f 



8748 



„ 2 



(5-1) 



where symm denotes (|2.24|) and ip(z) is the derivative of the logarithm of the Euler T-function. 
Other various quantities are defined by 



n 



(2) 
31 



, -1/(2) 



Sn(^) 



1 



V~3 



73) 



(5.2) 
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where Li n (z) is the polylogarithm function. We retain the notation of [13] here in defining the 

(2) (2) 

quantity S which is a linear combination of two harmonic polylogarithms, T~L\i and ■ More 
background to their appearance in the basic master integrals can be found in Appendix A of 
|13j . The theory for such harmonic polylogarithms was developed in |50| . These polylogarithms 
together with the other quantities such as £(3) and ln(3) emerge from the one and two loop 
master integrals, [Ml 071 0HJ 09] . Recently, the mathematics of these masters has been studied 
in the context of cyclotomic polynomials and harmonic polylogarithms in [51]. There an insight 
has been given for which particular polylogarithms and other such numbers will arise in the 
higher loop order master integrals. We also note that our convention in (|5.ip is that when a 
Lorentz index is contracted with one of the external momenta p or q then the Lorentz index is 
replaced by the momentum to compactify notation. 

We only include the expression for the Green's function containing 0\ since that for O2 
can be readily deduced by multiplying (I5.ip by 7 s (g) 7 s . Then the sum and difference of these 
two cases will give respectively the finite parts of the Green's function of the eigen-operators of 
7y(a). However, we have checked that both are in agreement in four dimensions. This is because 
while we have derived the finite renormalization required to restore anti-commutativity of 7 s in 
four dimensions that was essentially based on the two structures I <g> I and 7 s ® 7 s . Aside from 
the fact that the 3-quark operators are renormalizable, the reason for this is that these channels 
ordinarily contain the divergences in e. Therefore in keeping with the Larin method, [31J, they 
are used to define the finite renormalization matrix. In choosing this procedure it transpires 
that the coefficients in the other channels are in agreement when the conventions on the charge 
conjugation matrix, C, are respected. Next in extracting (15, ip from our d-dimensional expression 
in addition to (|2.15p we have to be careful in taking the four dimensional limit of generalized 
T-matrices where there are contractions with either or both of the external momenta p and q. 
So in addition to (|2.15p we have used 



1 (4) ^ 1 (4) qiivo 



= 6pqj ® 7 



1 (4) 59 1 (4) pqfw 



5 ® 7 5 

2„2„,5^„,5 n / \2„5 ^ „,5 



= 2pV7 ®7 - 2(pg)V®7 (5-3) 

=4 

for ()2.22p . The explicit values, (|2.24j) and (|2.25p . can be substituted in these general expressions. 
We also should comment on the structure of (I5.ip in the various channels. One approach to 
extract the finite part might have been to project the coefficients of each Lorentz structure. 
While we did not follow that line one can observe from the final result that there are seven such 
channels at two loops in four dimensions. In d-dimensions, prior to taking the limit to four 
dimensions after renormalization, there are more than seven channels. Moreover, it seems that 
not all possible structures are generated and it would appear that that strategy would require 
more effort than is necessary. Also internally it would require the manipulation of long strings 
of 7-matrices after applying the transpose of a set of these matrices due to (]2.2p . This would 
slow the symbolic manipulation programmes. 

While ()5.ip represents the evaluation of the Green's function to two loops in MS analytically, 
for practical purposes the numerical value is more useful. Therefore, we have evaluated (15. ip to 
seven decimal places and find 

(Vv(p)# O?)^ (-P - q)Oi 5(0)) = [1 + [0.9894261a + 0.9894261] a 

symm 

+ [41.5310566 + 6.7082190a + 2.8995053a 2 

- 3.9141771JV/] a 2 ] I a5 <g> 1^ 
- 1.6908864a 2 7^ ® 7 J 7 
+ [[0.5208683a + 0.5208683] a 
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[16.3956216 + 6.1021902a + 1.8664447a 2 

2] I a s ® 



1.2732336iV>] a . 



0.3372607a 2 



M 2 



+ [[0.2604341a + 0.2604341] a 

+ [9.4012326 + 3.0510951a + 0.9332223a 2 



0.5498054iV) ] a 2 ] — - 



+ [[0.5208683a + 0.5208683] a 

+ [17.4304130 + 6.1021902a + 1.8664447a 2 

- 1.1864222iV f ]a 2 ] 

+ [1.3720521 + 0.0868114iV f ] a 2 ^ — 

+ 0(a 3 ) (5.4) 



where we have used 



C(3) = 1.20205690 , E = 6.34517334 , $ f = 10.09559713 , 

ifj>" = 488.1838167 , s 2 (|) = 0.32225882 , s 2 {^j = 0.22459602 , 

s 3 (-) = 0.32948320 , s 3 (-) = 0.19259341 (5.5) 



27 V6 

as the input values in this exercise. Finally, we have not included the finite part of the Green's 
function (I2.26P to two or three loops since lattice measurements would require a zero momentum 
quark. This is a very difficult task numerically on the lattice. Moreover, only the two loop results 
of (12.22p are in an arbitrary linear covariant gauge since we restricted the three loop Mincer 
calculation to the Feynman gauge. 



6 General operator 

We now consider a generalization of the basic spin \ operators we have focused on so far. 
Recently, Krankl and Manashov, [32], have introduced the operator 

Of M = J JK 1%1?! J tf* (6-1) 

which has no contractions over the spinor indices and i, j and k are flavour indices. By con- 
sidering the renormalization of this basic operator and its mixing into operators of the same 
dimension they managed to derive two loop expressions for several other operator aside from 
the two considered in the previous sections. However, the anomalous dimensions recorded in [32] 
are not in the MS scheme. Despite this we have extended the results of [32] to three loops. This 
is straightforward as the symbolic manipulation programmes used to derive the MS renormaliza- 
tion of the (|, 0) operators was sufficiently general and hence adaptable to (16. ID . First, we recall 
the notation and formalism for the renormalization of (|6.ip . [32] . The bare operator mixes into 
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an infinite set of related operators which involve the generalized 7-matrices, T^\"^ n . Though at 
each order in perturbation theory the number of generated operators is finite. Specifically, 

n ijk _ 7 a'P'Y n ijk re. fj\ 

where the renormalization constant matrix Z™a ^ ^ is given by 

Z aa' /3 /3' -y-y' ^aa' 

(6.3) 

k 

Here the poles in e are contained within the function a mnp (e) where h = mnp is a label which 
indicates the basic Tr n \ -matrix structure and there is no sum over individual m, n and p but 
over the corporate label k, and 

00 (n.) 

a) = E^F • (6-4) 

n=l 

As there are three open spinor indices in f)6. 1 1) the T^-matrix structure is of the form, [32 1, 

G m np(ao'|^'|77') = F(m)aa' ® F („)££' ® . (6.5) 

We will omit spinor indices from this point and use tensor product notation as it is clearer. We 
have not included the Lorentz indices here but the contractions are across different T-matrices 
due to the antisymmetric property and there are no free Lorentz indices. It turns out that from 
explicit calculations the r( n )-structures appear in a symmetric form. To three loops, using the 
same notation as [32], these are 

Co = TooO , C2 = T220 + T202 + To22 , C4 = T440 + 1^404 + To44 

^6 = T660 + T606 + To66 , C222 = ^22 , C42 = 1^22 + 1^42 + 1^224 

C442 = T442 + T424 + T244 ; C444 = ^44 

C642 = T642 + 1^624 + 1^462 + 1^26 + 1^46 + 1^264 (6.6) 

where Cj retain the same spinor index structure as G. We have included C222 and C442 in this list 
as they appear at intermediate parts of the renormalization but are absent in the final expression 
in keeping with the expectation that the total number of 7-matrices should be divisible by four. 
To avoid any confusion the explicit contraction of the Lorentz indices in each definition of (16.6H 
is 
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(6.7) 



with the obvious permutation of m, n and p to define the forms in (|6.6p not listed in (|6.7p . In 
(16. 6h Cq, C2, C4 and C42 arise at two loops, [32], and Cg, €222, C442, C444 and C642 only appear 
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at three loops. That there are no other structures to this order is elementary to deduce from the 
fact that beginning with (|6.ip there are 4, 8 and 12 possible 7-matrices in each of the respective 
one, two and three loop Feynman diagrams. With the absence of free Lorentz indices and the 
antisymmetry (j6.6|) are all that survive. 

In |32] in order to ease the derivation of the anomalous dimension of ()6. 1 j) at two loops in 
four dimensions a relation was derived for C42 in d-dimensions which was 



C 



12 



3d(d - 1)C - 2(d - 3)C 2 



7^4 + -< 



(6.8) 



where the product of the Q is regarded as the multiplication of the constituent 7-matrices. 
While this is a relation in d-dimensions ultimately we will require the anomalous dimension in 
four dimensions and as noted in 1 32 1 then 
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In addition we have the similar but more trivial relations 
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(6.10) 

At three loops two new structures emerge, Cq and C64 2 , which involve the evanescent r^'"^ 6 
matrix. However, similar to (j6.8j) one can deduce that in d-dimensions 

C fi = - 12d(d- l)(2d- l)C n - 3(d - l)(7d- 24)C 2 - 6(2d - 5)C 4 

3C444 

: 642 = i2d(d - l)(2d - 7)C + 9(d 2 - 9d + 16)C 2 + 2(2d - 5)C 4 

3C 444 • (6.11) 



- 12d(d- l)(2d- 1)C 
+ 2(3d - 4)C 2 2 - ici + ^C 2 C 4 



2(3d-io)c 2 ! + ic| - ^C 2 C 4 



So these evanescent combinations can be expressed in terms of Cj which do not involve any 

p/ii.../i n w jt]j n > 5 p or the restriction to four dimensions we have the additional relation 

(n) - 



C444 
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(6.12) 



This follows trivially from the antisymmetry property. Given the presence of P 



possibility for this could have been 7 5 
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one 
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7 . It is easy to see that this is excluded when 



one examines the pattern of Lorentz indices in four dimensions. The use of the relations ()6.1ip 
can be viewed within the approach of [32] as a variation of the formalism of |30| which was 
introduced to include the effect of evanescent operators in the renormalization group functions. 



Equipped with these identities we have extracted the three loop anomalous dimensions (|6.ip 
in the MS scheme from the same three loop Mincer computation as in previous sections. 
Though in this case we do not contract the free spinor indices to produce a spin | operator. 
The full result is 
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(6.13) 



The two loop part is in exact agreement with [32]. Moreover, our two loop computation was 
carried out in an arbitrary linear covariant gauge and we observed the cancellation of the gauge 
parameter which provides an additional check. The three loop diagrams were computed in the 
Feynman gauge and the double and triple poles in e in the three loop renormalization constant 
satisfy the underlying renormalization group formalism as otherwise a finite expression would 
not have emerged. Therefore we are confident that (|6.13p correctly extends the result of |32j . 
If one restricted to four dimensions then the final term involving C444 would be absent and the 
corresponding four dimensional expressions for the remaining used. 
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Table 1. Values for the evaluation of the general anomalous dimension for various nucleons. 



Equipped with the general anomalous dimension we can extend the two loop results in the 
renormalization scheme of [32j for eigen-operators with specific spins and chirality. As indicated 
in [32] these relate to various operators in the literature. If we denote the label which the 
irreducible representations of the Lorentz group with two spins j and j by have then the 
eigen-operators are 0w>i). They can be written in the forms, [23, 
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where A 2 = 0. Here we have denoted right and left handed quarks by tpR = i(l + 7 5 )^ an d 
ip L = i(l — 7 5 )i/i respectively. To determine the anomalous dimensions of each operator from 
the general anomalous dimension we replace these four dimensional tensor product matrices Q 
by their eigenvalue under the Lorentz symmetry, |32j . These have been given in [32] but are 
summarized in Table 1. Hence we have 
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To assist with a comparison the numerical values are 
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+ [2.07407417V/ 2 - 47.94514927V/ + 618.3440216] a 3 + <3(a 4 ) 
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(6.16) 



Clearly the coefficients of the anomalous dimension derived from the generalized operator ap- 
proach of [32j are larger in value than our MS direct calculation. Therefore it would appear that 
the latter anomalous dimensions have a slower rate of convergence. 

Comparing the expressions for the (i, 0) pair of operators to our MS expressions we see that 
the one loop terms are the same. This is expected since that part of an anomalous dimension is 
scheme independent. The two and three loop terms are not the same. This difference is due to 
renormalization scheme dependence. While the subtraction method used is in principle the same 
in both cases since only the poles in e are removed into renormalization constants, it is in the 
derivation of these poles in the d-dimensional calculations where the differences arise. Moreover, 
in [32] the issue of having to handle 7 s in dimensional regularization is circumvented unlike 
our extension of the two loop MS computation of [20] where we completely reproduced that 
result. One issue relating to this concerns whether it is possible to derive (|4.8p from the general 
operator anomalous dimension, (16. 13H . A clue resides in the comparison of the expression for 
both chiralities of the spin 4 operators. If we compute the ratio of the anomalous dimensions of 
the chiralities for both spin | cases we find 
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As in [32] the first two terms are in agreement but differ now in the 0(a 2 ) terms. Moreover 
coefficient of the one loop /3-function appears in the discrepancy between the MS result of 
and (|6.15p . This can be explained by recalling that in the mapping of the general result C2 is 
replaced by 12. However, this is the four dimensional evaluation of did — 1) which is derived 
from the product r^r^)^- Expressing this in terms of e gives 
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Such a factor emerges from C2 when one projects (|6.13p onto the proton operator by using 
the formal contraction / <g> I. Therefore, including the contribution from the O(e) term within 
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the mapping used in [32] reproduces the corresponding two loop terms of (|4.8p . To extend 
this to the next order is certainly highly non-trivial. This is partly because there are more 
terms corresponding to additional operators in (I6,13|) but also due to the now hidden general 
evanescent operators as well as the 7 s issue. The operators are hidden in the sense that their 
effect in d-dimensions cannot be restored from the explicit expression in (I6.13p . More crucially, 
though, in order to proceed along these lines one is in effect repeating the actual calculation 
anyway which was used to derive (|4.8p in the first place. Indeed an analogous analysis for 4-fermi 
operators, [52], only serves to illustrate the large complexity of such a problem which is beyond 
the scope of the present article. 

7 Discussion. 

By way of concluding remarks we note that first we have extended the two loop MS renormal- 
ization of the 3-quark proton operator, [2D], to three loops. This is a technically more involved 
computation than [20] since the 7 s problem in dimensional regularization can no longer be 
treated passively at three loops. To accommodate this we have extended Larin's method for 
automatic symbolic manipulation programmes to operators which mix under renormalization. 
Similar features to [31] emerge in that the finite renormalization constant which is required 
to restore anti-commutativity in four dimensions is independent of the gauge parameter. We 
have indicated that this is true to all orders if the naive anomalous dimensions of the operators 
are independent of the gauge parameter which is the case for MS. By contrast in other mass 
dependent renormalization schemes this finite renormalization would be gauge dependent. The 
reasoning for this is that when a gauge invariant operator is renormalized in a mass dependent 
scheme its anomalous dimension depends on the gauge parameter. In addition what is apparent 
from comparing the various finite renormalization constants used to restore chirality in four 
dimensions both here and in [31] is that there is no universal finite renormalization. In other 
words one cannot merely extract a Z 5 from [31] and use it within another computation where 
the seed operator is not even present. While it may appear to be satisfactory at a particular low 
loop order for a Green's function it will lead to inconsistencies at higher loop order. In other 
words for each appearance of 7 s within an automatic symbolic manipulation computation one 
has to systematically treat 7 s in an adaptation of the Larin method. 

Moving away from the above general remarks concerning 7 s we note that the 3-quark operator 
we concentrated on was that relating to the proton or (^,0) in Lorentz spin notation. Ideally 
other spin operators are of interest and we have provided a first step in that direction by 
extending the recent analysis of [32] to three loops. This used a generalized operator approach 
which resolved the evanescent and 7 s issues from another angle. Though as we have remarked 
it is clear that the results are not in the MS scheme if one regards the earlier work of |31] as 
the true MS situation which we are in agreement with. To extend the 3-quark operator MS 
renormalization to other spin operators is not straightforward. This is because there is mixing 
between operators which is dependent on flavour symmetry. Our computational setup was 
designed purely for the proton case and will need to be extended to accommodate these other 
operators which is a topic we hope to return to later. The explicit definition of these operators 
and the relation to the flavour structure is given in [24. 39 . Whilst the motivation for this 
work is in relation to providing the perturbative structure of Green's functions to assist lattice 
matching in the high energy limit, in order to refine the understanding of proton structure will 
require an extension of our analysis in another direction. In essence this involves the treatment 
of operators with higher moments which manifest themselves in the decoration of (12, ID with 
covariant derivatives. For instance, a three loop renormalization of the first moment will require 
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the order of an additional 700 Feynman diagrams to be calculated. Again we hope to return to 
this in a later analysis. 
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